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ABSTRACT 

Physical mechanism responsible for high viscosity in accretion disks is still under debate. 
Parameterization of the viscous stress as aP proved to be a successful representation of this 
mechanism in the outer parts of the disk, explaining the dwarf novae and X-ray novae outbursts 
as due to ionization instability. We show that this parameterization can be also adopted in the 
innermost part of the disk where the adoption of the a-viscosity law implies the presence of 
the instability in the radiation pressure dominated region. We study the time evolution of such 
disks. We show that the time-dependent behavior of GRS 1915+105 can be well reproduced if 
a-viscosity disk model is calculated accurately (with proper numerical coefficients in vertically 
averaged equations and with advection included), and if the model is supplemented with (i) 
moderate corona dissipating 50% of energy (ii) jet carrying luminosity-dependent fraction of 
energy. These necessary modifications in the form of the presence of a corona and a jet are well 
justified observationally. The model predicts outbursts at luminosity larger than QAQMEdd, as 
required, correct outburst timescales and amplitudes, including the effect of increasing outburst 
timescale with mean luminosity. This result strongly suggests that the a-viscosity law is a 
good description of the actual mechanism responsible for angular momentum transfer also in the 
innermost, radiation pressure dominated part of the disk around a black hole. 

Subject headings: accretion, accretion disks — binaries: close — black hole physics — instabilities — 
stars: individual (GRS 1915+105) 



1. INTRODUCTION 

A standard, geometrically thin accretion disk 
(Shakura & Sunyaev, 1973) is known to be un- 
stable in the innermost, radiation pressure domi- 
nated regions (Lightman & Eardley 1974, Pringle, 
Rees & Pacholczyk 1974)). The presence of this 
instability is due to an assumption that the viscous 
torque is proportional to the total (gas+radiation) 
pressure. Since our understanding of the nature 
of viscosity in a disk is poor (for a review, see Pa- 
paloizou & Lin 1995) we can explore the viscosity 



mechanism indirectly, through modeling the con- 
sequences of the radiation pressure instability and 
comparing them to the observed variability of X- 
ray sources. 

On the theoretical side, the basic description of 
stationary solutions was completed by Abramow- 
icz et al. (1988) who found that radial advec- 
tion stabilize the disk at high accretion rates so 
a time dependent behavior was expected only for 
intermediate accretion rates. The disk evolution 
was first computed by Honma, Matsuoto & Kato 
(1991), and Szuszkiewicz & Miller (1998) pre- 
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sented computations of several consecutive out- 
bursts, thus confirming a limit cycle behavior of 
the disk. 

On the observational side, a perfect candidate 
has been found to study the details of radiation 
pressure instability. This is the X-ray source GRS 
1915+105 discovered by GRANAT observatory 
(Castro-Tirado, Brandt & Lundt 1992) and ex- 
tensively observed by the Rossi X-Ray Timing Ex- 
plorer (Belloni ct al. 1997a, 1997b). The source 
exhibits a superluminal jet (Mirabel & Rodriguez 
1994). During outbursts, its X-ray luminosity is 
dominated by a disk-like component (Belloni et 
al 1997a), and therefore the observed variability 
in this state should reflect the time behavior of 
the disk. The variations are well modeled phe- 
nomenologically by periodic changes of the disc 
inner radius (e.g. Feroci et al. 1999). The sce- 
nario, in which the radiation pressure instability 
is actually responsible for these variations was out- 
lined by Belloni et al. (1997c). The duration of the 
hard spectral state episodes in the disk-dominated 
epochs was recently identified with the viscous 
timescale of a standard, radiation pressure dom- 
inated disk (Trudolyubov, Churazov & Gilfanov 

1999) . 

However, the time evolution computations of 
Honma et al. (1991) and Szuszkicwicz & Miller 
(1998) based on radiation pressure instability 
strongly overpredict the amplitude of the outburst 
of GRS 1915+105. Therefore, Nayakshin, Rapa- 
port & Melia (2000) suggested that the standard 
viscosity law does not apply and they successfully 
reproduced the overall shape of GRS 1915+105 
RXTE lightcurve under an assumption of viscos- 
ity law of Taam, Chen & Swank (1997). In order 
to reproduce the observed minimum luminosity 
required for outburst behavior they also had to 
assume that the disk dissipates only 10% of the 
energy, with the remaining energy dissipated in 
a hot corona, which is not supported by observa- 
tions. 

In this paper we apply the standard viscosity 
assumption and nevertheless we successfully re- 
produce the RXTE ligthcurvc of GRS 1915+105. 
The key elements of our model, in comparison with 
Honma et al. (1991) and Szuszkicwicz & Miller 
(1998), are: (i) the presence of the outflow in 
a form of a jet (modeled after Nayakshin et al. 

2000) , (ii) the use of the appropriate coefficients 



in vertically averaged equations determined from 
the disk vertical structure, (iii) the presence of a 
moderately strong corona, dissipating 50% of the 
energy. 

2. METHOD 
2.1. Assumptions 

We study the time-dependent structure of a Ke- 
plerian vertically averaged accretion disk. 

We use a standard viscosity prescription, i.e. 
we assume that the viscous stress tensor is pro- 
portional to the total pressure, P (r rv = —aP). 

The angular momentum distribution of disk 
material is approximated as Keplerian so we ne- 
glect the problem of the transonic character of the 
flow close to the marginally stable orbit studied 
by Honma et al. (1991) and Szuszkicwicz & Miller 
(1998). Instead, we pay much more attention to 
various disk cooling mechanisms. 

We assume that the heat generated within the 
disk at any radius is either stored temporarily 
within this radius or is removed by: (i) radiation 
(ii) advection (iii) jet. Because the spectral obser- 
vations show the presence of a hard X-ray tail, we 
also assume the existence of a hot corona. 

The radiative cooling of the disk is calculated 
assuming that the disk is optically thick and radi- 
ates as a black body. The fraction of energy car- 
ried by advection is determined by standard equa- 
tions (e.g. Paczyhski & Bisnovatyi-Kogan 1981, 
Muchotrzeb & Paczyhski 1982, Abramowicz et al. 
1988). 

Fraction of energy, fj et , carried locally by a jet 
is parameterized after Nayakshin et al. (2000): 

Here the constant A is a model parameter describ- 
ing the strength of the jet and m is the local accre- 
tion rate at a given moment and radius, measured 
in units of the Eddington accretion rate: 

,v L Edd AirGMm H 



The assumed efficiency of accretion is i] = 1/16, 
as it results from the pseudo-Newtonian approxi- 
mation to disk accretion. 
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In our model this cooling mechanism is included 
in the energy balance, in opposite to the model of 
Nayakshin et al. (2000) where the jet was only 
used as an energy channel carrying a fraction of 
energy dissipated in the corona. 

The fraction of the energy dissipated in the 
corona is a free parameter of the model and in the 
present model we assume it takes a constant value, 
fcor = 0.5, independent from time and radius. 

The time evolution of the disk is governed by 
two equations: energy balance given in the follow- 
ing form 



dS 

£T — F n pn F rar j F n rii, Fqp± F r , 



(3) 



and the standard equation describing the time evo- 
lution of the disk surface density 

where the cooling by jet and corona is defined as 
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(6) 



Determination of algebraic relation between the 
quantities in Eq. 3 and 4 results from the con- 
tinuity equation, hydrostatic equilibrium and ra- 
diative transfer. It involves the determination of 
dimcnsionlcss coefficients which results from re- 
placement of the disk vertical structure with ver- 
tically averaged (or equatorial) quantities (e.g. 
Abramowicz et al. 1988, Honma et al. 1991). 

We determine those coefficients from the study 
of the vertical structure of a stationary disk at 
10 r g (see e.g. Dorer et al. 1996). Apart from 
standard ingredients, our stationary disk model 
included the effect of radial advection, energy 
transport by convection and appropriate bound- 
free opacities. The numerical code used for the 
computations was developed from the version of 
Pojmahski (1986) and subsequently modified by 
Rozariska et al. (1999). The coefficients, as 
defined by Muchotrzeb & Paczyfiski (1982) and 
Abramowicz et al. (1988) are: Bl = 0.8, S3 = 5.0 
and BA = 6.0, while the above authors used 
Bl = 0.67, B'i = BA = 6.0 and Honma et al. 
(1991) used Bl = 1.0, S3 = 16.0 and BA = 8.0. 



The importance of the model ingredients is il- 
lustrated by the stability curve of a stationary 
model calculated at 10 Rschw for 10M Q black hole 
mass and viscosity parameter a equal to 0.01 (see 
Figure 1). 

The triangles in Fig. 1 are obtained from the 
solution of the disk vertical structure, and dashed 
line is the stability curve obtained from the ver- 
tically averaged disk equations (with appropriate 
coefficients). In these two cases the disk corona 
and jet were neglected (A = and f cor = 0). The 
dotted line shows a modification of the stability 
curve by the corona (f cor = 0.5). The solid line 
shows the effect of the jet (A = 0.05). 

On the same plot (circles) we show the stability 
curve calculated similarly to the method adopted 
by Nayakshin et al. (2000), i.e. we computed 
the disk surface density from the vertically aver- 
aged disk structure as specified in that paper, but 
we assumed the standard viscosity law. For the 
Nayakshin et al. model the instability occurs for 
accretion rates as low as to ~ 0.03. Therefore, in 
order to stabilize the disk at low to the authors 
had to postulate a strong corona. In our model 
we expect the instability to appear for accretion 
rates higher than m ~ 0.1 at that radius and even 
higher than to ~ 0.2 if half of the gravitational 
energy is dissipated in the corona (f cor = 0.5). 

2.2. Time Evolution and Model Parame- 
ters 

We use the time evolution code originally devel- 
oped and described in detail by Smak (1984) for 
the study of cataclysmic variables, and modified 
by Sicmiginowska, Czerny & Kostyunin (1996) in 
the context of ionization instability in AGN (see 
also Hameury et al. 1998). 

We adopt the mass supply, to d , to the inner part 
of the disk as a model parameter and follow the 
time evolution of the disk under the radiation pres- 
sure instability in thermal and viscous timescale, 
as studied by Szuszkiewicz & Miller (1998). The 
complete model is given by the mass supply rate, 
to , viscosity parameter, a, jet efficiency factor, 
A (relating the fraction of energy carried by the 
jet to the local temporary accretion rate), and the 
(fixed) fraction of energy dissipated in the corona, 

fcor- 
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3. RESULTS 

We calculate the time evolution of the disk 
for a set of parameters adequate to explain the 
timescales and the amplitudes of the typical out- 
bursts observed in GRS 1915+105. In the follow- 
ing we assume that the black hole mass is equal to 
M = 10M© (estimated mass of the microquasar 
ranges from 7 to 33 M Q ), viscosity parameter 
a = 0.01 (as used by Nayakshin ct al. 2000 on 
standard branch), f cor = 0.5 (moderate corona) 
and A = 0.05 (mildly strong jet efficiency). 

The instability can occur only above a certain 
accretion rate limit for a given mass and viscosity. 
For our choice of parameters, the disk is stable 
when the accretion rate m is smaller than 0.16. 
For accretion rates above this limit the disk ex- 
hibits strong and regular outbursts. 

The amplitude of the global outburst depends 
mostly on the shape of the stability curve in the 
innermost part of the disk, which constrains the 
jet efficiency A. The duration of the limit cycle, 
on the other hand, is basically determined by the 
viscous timescale at the outer radius of the insta- 
bility zone in a stationary model, so it scales with 
the choice of the viscosity parameter a. 

In Figure 2 we show four lightcurves calcu- 
lated for four values of the external accretion rate. 
Other parameters were kept the same. The plot- 
ted lightcurves represent the bolometric luminos- 
ity of the disk. However, they may be to some ex- 
tent directly compared to the observed lightcurves 
expressed in counts/sec since the count rate in 
RXTE is strongly dominated by the soft energy 
band, determined mostly by the disk luminosity. 
The outbursts are almost periodic but their overall 
shape is significantly influenced by the mass sup- 
ply rate. When m is only slightly larger than the 
minimum value required for instability to operate, 
the outbursts last only ~ 100 s and the separation 
between them is relatively long. With an increase 
of m the duration of the bursts increases up to 
~ 1000 s although the amplitude does not change 
strongly. Similar results were obtained by Nayak- 
shin et al. (2000). Although their approach to 
the description of the disk structure was different 
from ours, it led to a similar stability curve and 
therefore outburst properties. 

By comparing our model to the GRS 1915+105 
behavior we conclude that the observed shape 



of the microquasar's lightcurves can be obtained 
by varying the mass supply rate. The observed 
outburst's amplitude constrains the jet efficiency. 
The energy carried away in the jet varies during 
the cycle: it is negligible at low luminosity but can 
reach ~ 15% — 20% during the outburst, for as- 
sumed A=0.05. The jet losses thus provide more 
efficient mechanism for stabilizing the disk than 
the advection and ensure that the luminosity of 
the source during outburst does not exceed the 
Eddington value. 

The viscosity parameter a does not influence 
the shape of the stability curve, therefore the am- 
plitude of the outburst remains the same as long 
as a is the same on the upper and lower branches. 
The change of a in this case results in the hor- 
izontal shift of the entire M — S curve. How- 
ever, this strongly affects the surface density and 
viscous timescale, so for small a the surface den- 
sity is higher and the viscous timescale is longer. 
The adopted value well represents the observed 
timescales in the microquasar. 

4. DISCUSSION 

We show that a physically viable instability due 
to radiation pressure in a standard disk model 
can well represent the time dependent behavior 
of GRS 1915 +105 if the advection and the pres- 
ence of a moderate corona (/ cor = 0.5) and a jet 
is included in the model. 

The mechanism operates for accretion rates 
above m ^0.16 which is in agreement with the 
observations of GRS 1915 +105 - the source does 
not exhibit any outbursts if the mean luminosity 
temporarily drops below 2.1 x 10 38 erg/s. The vis- 
cosity coefficient of order of 0.01 is appropriate to 
model the typical outburst duration as a function 
of the mass supply to the innermost part of ac- 
cretion disk. The presence of the jet is necessary 
to explain the observed amplitude of the outburst, 
since without the jet the lightcurve variations are 
too large (factor ~ 30). The role of advection and 
the jet in the presented parameterization are sim- 
ilar, but for the required value of the jet efficiency 
the jet losses dominate and the effect of advection 
is never strong. 

We conclude that the parameterization of the 
viscous stress as proportional to the total pressure 
well represents the true viscous stress properties 
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and the radiation pressure instability is a promis- 
ing model of the basic instability mechanism un- 
derlying the observed variability of this micro- 
quasar. The energy losses due to the jet are an 
essential ingredient of such model. Models based 
on modified viscosity law may allow for much lower 
jet efficiency. Observational constraints on the 
amount of energy carried by the jet will allow to 
distinguish between the two possibilities. 
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Fig. 1. — The stability curves calculated for M = 
1OM , a = 0.01 and R = lOR g . The triangles 
mark the solution with vertical structure model 
including advection and the circles mark the stan- 
dard analytic solution. The dashed line is the so- 
lution with scaling adopted to match the vertical 
structure model. The dotted line marks the so- 
lution obtained from our analytic approximation 
with the corona included (f CO r — 0.5), and the 
solid line shows the solution with the jet (A=0.05). 




Fig. 2. — The lightcurves calculated for M = 
10M Q , a — 0.01. and four values of the external 
accretion rate: 0.6 x 10~ 7 , 0.8 x 10~ 7 , 1.0 x 10~ 7 
and 2.0 x 10 _7 [M Q /j/r] . The jet ejection de- 
scribed by Equation (1) is included with coefficient 
A = 0.05 and the corona contribution is fixed at 
fcor = 0.5. 
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